Evolution of a Metal to Insulator Transition in Ca^zNa^CuC^C^, as seen by ARPES 



F. Ronning 1 *, T. Sasagawa 2 , Y. Kohsaka 2 , K.M. Shen 1 , A. Damascelli lt , C. Kim 1 - 5 *, T. Yoshida 1 , N.P. Armitage 1§ , 

D.H. Lu 1 , D.L. Feng x t, L.L. Miller 3 , H. Takagi 2 , Z.-X. Shen 1 
1 Department of Physics, Applied Physics and Stanford Synchrotron Radiation Laboratory, 
Stanford University, Stanford, CA 94305, USA 
2 Department of Advanced Materials Science, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
3 Department of Physics, Iowa State University, Ames Iowa, 50011 
(February 1, 2008) 

We present angle resolved photoemission (ARPES) data on Na-doped Ca2Cu02Cl2- We demonstrate 
that the chemical potential shifts upon doping the system across the insulator to metal transition. 
The resulting low energy spectra reveal a gap structure which appears to deviate from the canonical 
d x 2_ y 2 oc \cos(k x a) ~ cos(k y a)\ form. To reconcile the measured gap structure with d-wave super- 
conductivity one can understand the data in terms of two gaps, a very small one contributing to 
the nodal region and a very large one dominating the anti-nodal region. The latter is a result of 
the electronic structure observed in the undoped antiferromagnetic insulator. Furthermore, the low 
energy electronic structure of the metallic sample contains a two component structure in the nodal 
direction, and a change in velocity of the dispersion in the nodal direction at roughly 50 meV. We 
discuss these results in connection with photoemission data on other cuprate systems. 



I. INTRODUCTION 

The close proximity of an antiferromagnetic Mott insu- 
lating phase suggests that this insulator may be the nat- 
ural starting point for understanding the unusually high 
superconducting transition temperatures as well as the 
unconventional normal state properties 1 in the cuprates. 
The insulator to metal transition in strongly correlated 
materials is a well studied fundamental problem in solid 
state physics in its own right. Even by limiting ourselves 
to a one band Hubbard model, which may be an appro- 
priate starting point for the half-filled cuprates, 2 there re- 
main several ways one could imagine the electronic struc- 
ture evolving upon doping (see figure 1). One scenario is 
the chemical potential shifting to the top of the valence 
band similar to the case of a band insulator. 3 An alter- 
native scenario is states being created inside the Mott 
gap upon doping. 4,5 These are the two most commonly 
discussed pictures with regard to the cuprates. 

ARPES studies of the Laa-zSr^CuC^ (LSCO) sys- 
tem show that spectral weight is transfered to low en- 
ergy states, which lie inside the Mott gap, upon hole 
doping. 4-7 In contrast, in the P^S^Cai-aY^C^Os+s 
(Bi2212) and Nd 2 _ 2; Ce 2; Cu04 (NCCO) systems there is 
work which shows that the chemical potential shifts in 
addition to a transfer of spectral weight. 8,9 A poten- 
tial problem for the Bi2212 system is the difficulty to 
obtain deeply underdoped samples. The oxychlorides, 
Sr 2 Cu02Cl2 and Ca2Cu02Cl2, are single layer cuprates 
that have provided the best photoemission data on the 
undoped insulating state to date. Thus the oxyhalides 
are ideal systems to study the metal to insulator transi- 
tion with photoemission, as high quality single crystals, 
which cleave as easily as the Bi2212 system, can now be 
grown through the metal to insulator transition. 
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FIG. 1. Schematic doping evolution of an insulator, a) 
Upper (UHB) and lower (LHB) Hubbard bands for a Mott 
Insulator, b) Result of doping holes into a band insulator, c) 
Result of creating states inside the gap with doping, d) A 
scenario combining the ideas of b) and c) . 

Recently, we have reported that the chemical potential 
shifts upon doping the insulator, Ca2Cu02Cl2(CCOC), 
with Na. 10 Here we present the full body of evidence 
which supports arguments for such a shift of the chemi- 
cal potential in the C^-zNazCuC^Cko system. We ex- 
amine the valence band and photon energy dependence 
of Na-doped CCOC and compare it with that found in 
the insulator. By making a detailed comparison of the 
low energy spectra we unambiguously show that the high 
energy pseudogap is in fact a result of the rf-wave-like 
dispersion seen in the insulator, as had been suggested 
previously. 11,12 

As we will show, residual effects from the insulator 
strongly distort the "band structure" in these deeply un- 
derdoped samples resulting in a gap structure along the 
LDA-predicted Fermi surface that strongly deviates from 
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the canonical d x 2_ y 2 form. In order to reconcile this re- 
sult with c?-wave superconductivity, we conclude that the 
measured gap structure is a consequence of probing two 
different regimes. One regime is concerned with a small 
unresolved gap along a hole pocket about (ir/2, n/2). The 
second regime is in the vicinity of (n, 0) where a large gap, 
which we identify as the pseudogap, is identified with the 
undoped insulator. This picture is consistent with the 
observed chemical potential shift. 

The evolution of the electronic structure, however, is 
not solely a simple shift of the chemical potential since 
the spectral lineshape also evolves with doping. We fur- 
ther investigate here the two component structure found 
in the nodal direction of Ca2- £C Na a; Cu02Cl2. 10 We show 
that this structure is not directly related to T c , and from 
a MDC (Momentum Distribution Curve) analysis we find 
a "kink" in the dispersion which has the same energy as 
the " dip" which divides the two components in the EDCs 
(Energy Distribution Curves) . The observation of a kink 
near 50 meV in this system supports the observed uni- 
versality of this type of behaviour in the cuprates. 13 Fi- 
nally, we discuss the similarities and differences between 
photoemission results from other cuprates and Na-doped 
CCOC. 



5 GPa gave a Na concentration of 10% and 12%, respec- 
tively. Crystals with a maximum dimension of 1.5 • 1.5 • 
0.1 mm 3 were obtained. The Na content was estimated 
by comparing the c-axis lattice constant with powdered 
ceramic samples of known content. 15 Magnetization mea- 
surements of T c gave results consistent with the assigned 
doping concentrations. The Na concentrations used in 
this study were x=0.10(T c =13) and x=0.12(T c =22 K). 
Since superconductivity does not occur until x^0.08 in 
the Na-doped CCOC system, 14 ' 15 we consider our sam- 
ples to lie in the heavily undcrdoped regime. As a refer- 
ence, the optimal T c in powder samples has been found 
to be 28 K. 15 ARPES measurements were performed on 
beamline 5-4 of SSRL. The samples were easily cleaved in 
situ and measured at the indicated temperatures in the 
figure captions. The angular resolution was 0.29°, while 
the energy resolution was always <30 meV. When exam- 
ining the detailed low energy excitations (ie after figure 
7) the energy resolution was improved to <15 meV. 

III. DATA AND RESULTS 
A. Valence Band 



II. EXPERIMENTAL 

Na-doped CCOC single crystals were grown by a flux 
method under high pressure, by using a cubic anvil type 
pressure apparatus. 14 A powder mixture of Ca2Cu02Cl2, 
NaCl, NaC10 4 , and CuO (1:0.2:0.2:0.1 molar ratio), 
sealed in a gold tube, was heated to 1250°C and then 
slowly cooled to 1050°C under high pressure. Uniform 
samples at various Na concentrations, x, were achieved by 
using different pressures during the synthesis. 4 GPa and 
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FIG. 2. a) ARPES data reveals three features, at roughly 
9, 10, and 12 eV, in Na-doped CCOC not present in the pure 
sample. The time elapsed between the initial cleave and the 
time the spectra were recorded is indicated, b) illustrates, 
with a different cleave, how the features vanish within 90 
minutes of the initial cleave. The photon energy was 25.5 
eV and the temperature was 20 K. The Fermi cutoff on the 
low energy spectral weight is too small to be visible in this 
intensity scale. 



We begin with the valence band spectra of Na-doped 
CCOC in figure 2, which illustrates a peculiarity of this 
system. Namely, from 8 to 13 eV binding energy there 
exists several pronounced features which rapidly vanish 
with time. The spectra from the cleave where these fea- 
tures were most prominent is shown in the top panel. 
Features are visible at 9, 10, and 12 eV which are not ob- 
served in the parent compound, CCOC. Figure 2b reveals 
that within 90 minutes following the initial cleave, the 
additional features from 8 to 11 eV have vanished. The 
feature at 12 eV also vanishes (not shown). For bind- 
ing energies less than 7 eV no change in the spectra was 
observed over these time scales where particularly close 
attention was paid to the near Ep features, less than 
1 eV binding energy The photon energy was adjusted 
to ensure the high binding energy features were not de- 
rived from a contribution of second order light from the 
monochromator. Comparing the spectra in panels a) and 
b) suggests a dispersive nature of the additional peaks 
not seen in the insulating compound. However, it was 
also observed that these features shift to lower binding 
energy as they age. We tentatively assign the origin of 
these features from 8 to 13 eV binding energy to flux 
inclusions, such as NaCl and NaC104 which are used in 
excess during the synthesis; 16 subsequently the surface of 
these flux inclusions are passivated by the photon follow- 
ing the cleave. Despite the uncertain nature of these high 
energy features, the fact that the low binding energy re- 
gion of the spectra have no variation over the short time 
scale in which the high energy features vanish, gives us 
confidence that our low energy results on the Na-doped 
CCOC crystals are intrinsic. 
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FIG. 3. Valence band spectra along the nodal direction for 
a) Ca 2 Cu0 2 Cl 2 and b) 10% Na-doped CCOC. £ 7 =25.5 eV, 
T=200 K and 20 K for x=0 and x=0.1, respectively. 

Figure 3 compares the valence band spectra along the 
nodal direction for x=0 and x=0.1. The data on the 
insulator are consistent with previous reports on a simi- 
lar material (S^CuO^C^), 17,18 and the data for the Na- 
doped Ca 2 Cu02Cl2 were reproduced on multiple cleaves. 
The position of the valence band in the insulator with re- 
spect to the chemical potential can vary as much as 1 eV 
which is believed to be a result of pinning the chemical 
potential at different impurity levels that exist for the 
different cleaved surfaces. It is important to note that 
the valence band data of the insulator shown in figures 3 
and 4 are at the minimum possible binding energy, since 
if they are shifted further towards the chemical potential, 
spectral weight from the tail of the valence band peaks 
would be pushed above the chemical potential falsely in- 
dicative of a metallic nature. The low energy excitations 
within 1 eV of the chemical potential, which the remain- 
der of this paper will focus on, are barely visible on this 
broad scale. Generally, the two data sets are quite sim- 
ilar, although differing relative intensities for a few of 
the bands about (0,0) can be seen. The other difference 
is that the insulator data are sharper than those of the 
metal. Perhaps this is an indication that the surface of 
Na-doped Ca2Cu02Ci2 is not as well ordered as the pure 
sample grown under atmospheric conditions. Figure 4a 
plots the second derivative of the above spectra, which 
can be used to more clearly identify the dispersion of 
various features. This method agrees with the dispersion 



FIG. 4. (color) Valence band comparison of Ca2Cu02Ci2 
and 10% Na-doped CCOC. a) is a second derivative plot of 
the EDCs presented in figure 3 retaining only points with 
negative curvature, and b) maps the dispersion of the vari- 
ous features as seen by eye. The measurement conditions are 
identical except for the temperature of the insulator which is 
raised to avoid electrostatic charging during the photoemis- 
sion process. _B 7 =25.5 eV, T=200 K and 20 K for x=0 and 
x=0.1, respectively. The data are consistent with a shift of 
the chemical potential with doping as discussed in the text. 

extracted by the peak positions of the various bands in 
the insulator and the metal as shown in panel 4b. The 
resolved bands appear to have shifted to lower binding 
energy with doping by roughly 300 meV, as would be ex- 
pected if this were a simple band material. One apparent 
exception exists from the band at 5.6 eV binding energy 
in the metallic sample near (0,0) in figure 4b. However, 
due to the broadness of the features in the metal, it is 
likely that this band is a superposition of several peaks, 
resulting in this apparent discrepancy. Note that in the 
insulator several bands are resolved in this same region. 
Thus, the overall picture from the valence band data is 
highly suggestive of a chemical potential shift with dop- 
ing. The most convincing demonstration for this is at 
(■7T, 7r) where the peaks for both concentrations are most 
clearly resolved. Here, the peak positions indicate that 
the insulator is shifted to higher energy by at least 300 
meV relative to the metal. The specific value will depend 
on which insulating data set is used, since the chemical 
potential in the insulator can vary substantially as men- 
tioned above. We again emphasize that the current bind- 
ing energy of the insulator could not be found at smaller 
binding energies without spectral weight being shifted 
unrealistically above the chemical potential. Thus, nei- 



ther charging of the insulator nor the variability of the 
chemical potential in the insulator can be responsible for 
the observed shift upon doping 

B. Near Ef Comparison 

To further support this idea, figure 5 presents spec- 
tra near Ep along the entire (0,0) to (n, 0) and (0,0) 
to (it, 7r) lines, where all the spectra from this insulat- 
ing data set have been shifted by 750 meV so that the 
centroid of the feature at (tt/2,tt/2) is at the chemical 
potential. The value of 750 meV is specific to this par- 
ticular data set, but is 300 meV more than allowed by 
the uncertainty in the position of the chemical potential. 
Along (0,0) to (n, 0) the spectra are remarkably similar. 
In many instances fine details of the line shape match ex- 
tremely well. For example, a second broad component in 
the electronic structure can be observed at roughly 600 
meV higher binding energy from the first feature creat- 
ing a heavily asymmetric lincshape. 19,20 This indicates 
that upon doping the chemical potential simply drops to 
intersect the top of the valence band similar to the case 
of an ordinary band material. This would naturally ex- 
plain the large high energy pseudogap seen at (it, 0) as 
a remnant property of the insulator as first suggested by 
Laughlin. 11 ' 12 

Along (0,0) to (ir, ir) the situation is less clear. As 
the feature moves towards lower binding energy a sharp 
Fermi cutoff appears in the metallic samples, as antici- 
pated by the large peak at the Fermi level from the ener- 
getically shifted spectra of the insulating sample. How- 
ever, the match between the spectra becomes increas- 
ingly worse as (it/2,it/2) is approached. Perhaps this 
simply indicates that a Fermi crossing has occurred be- 
fore (n/2, 7r/2), and thus the weight at (ir/2, ir/2) is sup- 
pressed relative to the insulator or that some spectral 
weight transfer in addition to the shift of the chemical 
potential has occured. To see if the suppression was due 
to matrix element effects we also changed the polarization 
to maximize the nodal direction cross section, but it re- 
mained weaker. This suppression possibly indicates that 
increased scattering which could effectively smear out the 
k resolution, and would have the greatest effect where 
the dispersion is the steepest (near (tt/2, it/2)), has oc- 
cured. However, this would require a nontrivial small an- 
gle scattering mechanism, which does not simply produce 
an angle independent background. Also, the cleavability 
of the Na-doped compounds is the same as the parent 
insulator, and laser reflections from the sample indicate 
flat surfaces, which argue against angle averaging as the 
cause for the observed suppression of weight. Typically, 
highly angle dependent valence band spectra as shown in 
figure 3 are indicative of a good surface, although the Na- 
doped valence band features are somewhat broader than 
in the case of the insulator. Pothuizcn et al. noted 
that in SCOC the lineshapc of the 2 eV binding energy 



peak at (ir, it) is identical to the lincshape of the Zhang- 
Rice singlet at (ir/2,ir/2). Given such a correlation a 
broadened low-energy spectral lineshape at (tt/2,tt/2) in 
Na-doped CCOC is perhaps expected when the valence 
band at 2 eV has been broadened, as observed here (see 
figure 3). This suggests that the small differences be- 
tween the doped system and the half filled insulator may 
indeed result from a difference in sample surface quality. 

Considering that the features of hole doped CCOC 
track the dispersion of the features in the half filled insu- 
lator, one might wonder exactly how similar are the elec- 
tronic states for the two doping levels. A good way to test 
this is by analyzing the photocmission matrix elements. 
Under the dipole approximation the photoemission inten- 
sity can be written as / oc f\A- p^i)] 2 where ^V/ are 
the initial and final state wavefunctions, and A ■ p is the 
perturbing Hamiltonian. 21 Practically, this means that 
the cross section of the photoemission process will depend 
on the photon energy and experimental geometry as well 
as the wavefunctions themselves. Thus, if two wavefunc- 
tions are to be considered as the same, then they should 

r -> (rc.a) r -> 
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FIG. 5. EDCs of Ca 2 Cu0 2 Cl 2 (dashed line) (from ref 12 ) 
are shifted by a constant in energy and compared with EDCs 
of 10% Na-doped CCOC (solid line). The overlap is extremely 
good with the exception of the features near (it/2,-k/2). The 
data are normalized at high binding energy for comparison. 
The slight rise in spectral weight above Ef observed in some 
EDCs of CCOC is due to the presence of a core level excited 
by second order light. This is not present in the metallic 
sample where the second order light contribution is heavily 
suppressed by the use of a normal incidence monochromator 
as opposed to the grazing incidence monochromator used for 
the insulating data. £ 7 =25.5 eV, T=100 K and 20 K for x=0 
and x=0.1, respectively. 
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FIG. 6. Photon energy dependence along the nodal direc- 
tion for Ca2-iNa^Cu02Cl2 for x=0.1 and x=0. For each 
doping, the same k range is shown, while the photon energy 
is indicated in each panel. T=20 K and 293 K for x=0.1 
and 0.0, respectively. Notice, that the modulation of inten- 
sity varies with photon energy in a similar fashion for both 
dopings. 



have the same dependence on the experimental condi- 
tions. We have found that CCOC experiences a fairly 
dramatic change in the modulation of intensity along the 
nodal direction when changing the photon energy from 
16.5 to 17.5 eV. In figure 6, we compare the EDCs along 
the nodal direction of a 10% Na-doped CCOC sample 
and an undoped CCOC sample using 16.5 eV and 17.5 
eV photons. Wc notice that the change in the modulat- 
ing intensity on going from i? 7 =16.5 to 17.5 eV is similar 
for both samples. Namely, for 16.5 eV photons there is 
much more intensity at (tt/2,tt/2) and higher k values, 
which vanishes on going to 17.5 eV photons. Even the 
existence of the second feature at 600 meV higher bind- 
ing energy from the lowest energy feature can be seen 
in both samples along the nodal direction. These re- 
sults show that the gross features of the wavefunction of 
the electronic states within 2 eV of the chemical poten- 
tial discussed thus far are remarkably similar. The pecu- 
liar photon energy dependence shown in figure 6 will be 
discussed separately with a much larger body of data. 22 
Along with the fairly rigid shift in dispersion of the low 
energy states and the valence band, the lincshapc com- 
parison and photon energy dependence show conclusively 
that the chemical potential indeed shifts upon doping the 
half filled insulator, Ca 2 Cu02Cl2, with Na. 



C. Metallic Nature of Na-doped Ca 2 Cu0 2 Cl2 

A naive expectation of a simple chemical potential shift 
upon doping is the creation of a small hole pocket cen- 
tered at (tt/2, it/2). However this contrasts with the ma- 
jority of ARPES results from the optimal and overdopcd 
regime which find a large Fermi surface centered about 
(tt, 7t) as predicted by band theory. 23 We shall thus take a 
closer look at the near Ef electronic structure to see if the 



simple picture of a chemical potential shift as illustrated 
in figure lb is correct. In other words, is the electronic 
structure unchanged with the exception of the position 
of the chemical potential? Figure 7 presents a sampling 
of the EDCs at selected k points. The cuts shown in (a) 
through (d) are taken parallel to the nodal direction. For 
several of the cuts close to the nodal direction a broad 
feature starts at high binding energy and approaches the 
Fermi energy with increasing k, where a sharp Fermi cut- 
off of the spectra is indicative of a Fermi surface crossing. 
We will return to the presence of a peak-dip-hump like 
structure near the nodal Fermi crossing later. The Fermi 
crossings, which have been determined by the spectra 
with the minimum leading edge midpoint (the binding 
energy where the spectral intensity is half of the intensity 
at the peak position), have been highlighted in the figure. 
We note that while a minimum value of the leading edge 
midpoint identifies the spectra which corresponds to &f, 
the value itself will be pushed back to higher binding en- 
ergy if this portion of the Fermi surface is gapped. As 
the cuts approach (tt, 0) the spectral features become less 
pronounced and it quickly becomes impossible to iden- 
tify a Fermi crossing as one can no longer clearly identify 
a peak position to meaningfully extract a leading edge. 
This is best illustrated in panel (e) where spectra along 
the ostensible Fermi surface predicted by band theory are 
presented. Near the node, there is a clear Fermi cutoff at 
E^, but this feature disappears quickly as one approaches 
the antinodal region. 

One can also visualize the data with an intensity map 
at Ef, obtained by integrating each EDC by a ±4 meV 
window about E^. The highest intensity contour, shown 
in figure 8a, forms what appears as a small arc. 24 The 
open circles indicate where a clear Fermi crossing can be 
observed as in the cuts shown in figure 7(a) through (c). 
Note this ends as the Fermi surface crosses the antifcr- 
romagnetic zone boundary. By assuming that the Fermi 
surface can not abruptly end in the middle of the Bril- 
louin zone, one is left with two possibilities. Either the 
Fermi surface still lies along the one predicted by band 
theory and has been heavily gapped as one approaches 
(tt, 0) or the Fermi surface is a small hole pocket cen- 
tered on or near (tt/2,tt/2). Wc investigate these two 
possibilities below. 

Let us first attempt to push the assumption that the 
large LDA-like Fermi surface holds for this heavily under- 
doped sample. We separate the large Fermi surface into 
two regions, as illustrated in figure 8b. In Region I the 
broad spectral function has a sharp Fermi cutoff, which 
clearly defines kp, while in region II one must extrap- 
olate the Fermi surface from the Fermi crossings found 
in region I assuming a large Fermi surface that does not 
terminate in the middle of the Brillouin zone. This is 
shown by the solid dots in figure 8. Note that in the 
region near (tt, 0), variations in the precise shape of the 
Fermi surface will not affect the following analysis as the 
spectra have little variation which is apparent from the 
EDCs in figure 7d. 
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FIG. 7. Selected EDCs of Na-doped Ca 2 Cu0 2 Cl 2 taken at 
the k points shown by solid lines in the cartoons. The lines in 
the Brillouin zone indicate all the k points where spectra were 
taken. The bold spectra in (a) through (c) indicate fci?- No 
spectra are highlighted in (d) as no Fermi crossing is observed. 
Spectra along the Fermi surface predicted by band theory are 
shown in (e). The dots in (d) and (e) indicate a high energy 



feature as discussed in the text. 
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FIG. 8. (color) (a) Ef map obtained by integrating the 
x=0.10 spectra by ±4 meV about Ef. The open dots in- 
dicate the Fermi crossing, while the solid dots extrapolate 
the Fermi surface to the Brillouin zone edge as would be ex- 
pected in band theory, (b) an illustration indicating the area 
of the resulting Fermi surface pocket obtained by symmeteriz- 
ing the observed Fermi crossings about the antiferromagnetic 
zone boundary (dashed gray line). By counting the number of 
electrons in the enclosed area a doping level is found which 
is indicated in the figure. Red reproduces the points in (a) 
for x=0.10, while blue are results from a sample with x=0.12. 
The data naturally fall into two regimes: I(II) distinguishes 
the region along the underlying LDA Fermi surface where a 
clear Fermi crossing does (does not) occur. 



Figure 9 plots the leading edge midpoint for spectra 
along the large Fermi surface just discussed relative to 
the value at 0=45° (the nodal position for a pure d x a_ y a 
functional form). In this way one can determine the gap 
structure along the Fermi surface. The open symbols in- 
dicate the leading edge midpoint for the Fermi crossing 
spectra found in region I. However, in region II the lack 
of a clear leading edge midpoint prevented us from iden- 
tifying a Fermi crossing in the first place. Hence, we need 
an alternative method to describe the dispersion. Near 
(it, 0) an energy scale (albeit a broad one) still exists in 
the spectra as indicated by a dot in figures 7 (d) and (e). 
Thus we can blindly normalize the maximum of the spec- 
tra below 300 meV to 1, and extract a leading edge mid- 
point via this simplified approach. This is what is shown 
by the filled symbols in figure 9. For a simple spectral 
function this procedure would give the correct gap value. 
In our case, one can see that close to the node this pro- 
cedure roughly matches our more careful analysis, but 
it also allows us to extend our gap analysis beyond the 
regime where a clear Fermi crossing has been identified, 
albeit in a crude way. The solid black line is the expec- 
tation for the canonical d-wave functional form. As men- 
tioned above, in the optimal and overdoped Bi2212 sam- 
ples a simple spectral function and a large Fermi surface 
is observed, and the d-wave functional form well describes 
the experimental gap structure. 25 However, in the heavily 
underdoped regime studied here, there is a strong devia- 
tion from the d-wave functional form. In order to recon- 
cile this gap structure with d-wave superconductivity 26 
our assumption of a large Fermi surface centered at (71", n) 
must break down. In this case, figure 9 is not physically 
significant as it maps the energy along two physically dif- 
ferent regions as defined in figure 8b, each with its own 
gap structure. In other words, by assuming a d-wave 
gap structure we have shown that only a fraction of the 
spectra used for figure 9 actually lie on the Fermi surface. 

Region I does indeed represent a portion of the Fermi 
surface about (jt/2, tt/2) which may contain a very small 
gap, while region II does not contain a Fermi crossing, 
and the resulting large gap identifies a larger energy scale. 
We have just ruled out the possibility that the Fermi 
surface in this sample is a large hole pocket centered 
at (7r, 7r). If we maintain that the Fermi surface can 
not abruptly end in the middle of the Brillouin zone we 
must conclude that the Fermi surface forms a small hole 
pocket, as naively expected from the observed chemical 
potential shift. The question then becomes, how does 
such a Fermi surface develop a gap consistent with d-wave 
superconductivity. Given the experimental resolution we 
can not comment on whether or not the flat region in 
figure 9 would itself fit a \cos(k x a) — cos(k y a)\ functional 
form, which is the simplest expectation for d-wave pair- 
ing, although models exist with small hole pockets and 
d-wave pairing for which this is not the case. 28 We note 
that for a superconductor with a transition temperature 
of 13 K, the maximum of the gap from BCS theory is 
expected to be less than 2 meV at the antinode and even 
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FIG. 9. The leading edge midpoint of the spectra along the 
Fermi surface illustrated in figure 8 relative to the value at 
the rf-wave nodal position. The open symbols are the result 
when a clear Fermi cutoff can be identified, while the solid 
symbols are obtained by normalizing the maximum of the 
low energy spectra to 1, as discussed in the text. The solid 
line is a plot of cos(2tp), the d-wave functional form. The large 
discrepancy from the canonical d-wave functional form when 
attempting to describe the data in terms of a band theory-like 
Fermi surface centered at (tt, tt) suggests that it is better to 
understand the data in terms of two regimes, labeled I and II 
in figure 8. 

smaller along the observed arc, which is too small for 
our experimental resolution and is not the objective of 
this study. Perhaps this gap could be measured with an 
improved spectrometer. 

By symmeterizing the hole pocket about the antifer- 
romagnetic zone boundary one finds an area of 5.7%, 
which implies a doping level of 11.4% due to the electron 
spin. This is in rough agreement with the hole doping of 
10% determined by comparing T c and lattice constants 
with results on powdered samples. 15 Again we note that 
the hole pocket is also not necessarily symmetric about 
the antiferromagnetic zone boundary, 28 which could sub- 
stantially modify these values. The idea of a small hole 
pocket about (7r/2, 7r/2) in Na-doped CCOC is also con- 
sistent with the observations of shadow bands 27 at higher 
binding energy. 10 Presumably, the coherence factors re- 
sponsible for the shadow band are very weak making it 
difficult to distinguish a second Fermi crossing from the 
background at ~Ep- We stress that while our findings are 
quite consistent with the small Fermi pocket picture, we 
did not see the shadow Fermi surface at zero energy even 
though we do observe a shadow band at higher energies. 

Meanwhile, the observation of a chemical potential 
shift in this material from the insulating state and the 
photon energy dependence of the lineshape emphasize 
that the features near (tt, 0) (region II) in the Na-doped 
Ca2Cu02Cl2 sample are a direct result of the magnetism 
which dominates the antiferromagnetic insulator. To il- 
lustrate that the energy scale near (tt, 0) is the same as 
that of the insulator, recall figure 5 where the spectra 
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FIG. 10. Comparison of EDCs at selected k points of 
10% Dy-doped Bi2212 (taken from ref 31 ) and 10% Na-doped 
CCOC. The Na-doped CCOC system suggests a much larger 
pseudogap than in Bi2212, both for the low energy pseudogap 
(thick bars) and the high energy pseudogap (triangles) . 

of the insulator matches that of the Na-doped sample at 
(tt.0). 

We have also repeated the above gap structure analysis 
on a sample with x=0.12 (T c =22 K), and some of the re- 
sults are presented in blue in figure 8 and as circles in fig- 
ure 9. The data reproduce the above analysis with a few 
small quantitative differences. The Fermi arc increases 
in size (shown in figure 8b), consistent with an increase 
in doping, while the binding energy near (tt, 0) reduces 
slightly which demonstrates that this energy scale is de- 
creasing with increasing doping. This is also expected 
as we have identified this energy scale with the mag- 
netism of the insulator whose strength should diminish 
together with the antiferromagnetic correlations as holes 
are added to the system. This has also been observed 
in previous studies on the high energy pseudogap. 29 We 
should note that the hidden message of describing the 
Fermi surface in the underdoped regime as a small pocket 
about (7r/2,7r/2) is that the system must at some point 
transform upon doping to the large Fermi surface cen- 
tered about (tt, tt) observed in the overdoped regime. In 
measuring the volume of the pocket for the x=0.12 sam- 
ple we find the doping level to be 21%. The fact that this 
analysis begins to show a deviation from a doping level of 
x perhaps begins to identify a transition to a state where 
the doping level should be described as 1-x (the expected 
volume of the large Fermi surface centered at (tt, 7r)). 30 

In comparison with other cuprates, the lack of a Fermi 
cutoff in the EDCs would normally be identified as an 
extremely large pseudogap. In an attempt to quantify 
this suppression of weight, we compare 10% Na-doped 
CCOC to 10% Dy-doped Bi2212 31 in figure 10. 10% Dy- 
doped Bi2212 is underdoped with a T c of 65 K giving it 
a T c /Max(T c ) of 0.72, while 10% Na-doped CCOC has a 
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ratio of 0.46. A comparison with 17.5% Dy-doped Bi2212 
(not shown) with a T c /Max(T c ) ratio of 0.28 yielded a 
similar comparison. One clearly observes the loss of spec- 
tral weight at the Fermi energy for the Fermi momentum 
of (it, 0.3n) when compared to the respective Fermi cross- 
ings along the nodal direction near (0.47T, OAw). There 
are two common methods for characterizing the pseu- 
dogap with ARPES. The first, is by the shift of the 
leading edge midpoint of a spectra to higher binding en- 
ergy relative to the chemical potential as indicated by 
the dashed lines. This is typically referred to as the low 
energy pseudogap. For Dy-doped Bi2212 this is reason- 
ably well defined, and gives a value of 30 meV. For the 
Na-doped CCOC sample, we find a value of roughly 50 
meV at (71", 0.37r). However, as previously noted, the lack 
of a Fermi cutoff near (it, 0) implies that the low-energy 
pseudogap is not well defined in the Na-doped CCOC 
spectra. As a result, we are inclined to characterize the 
pseudogap by what is known as the high energy pseudo- 
gap which identifies a larger energy scale in the spectra 
about (-7T, 0), and is sometimes referred to as a "hump" . 29 
Here, we find a high energy pseudogap of roughly 300 
meV for Na-doped CCOC, compared with 120 meV for 
Dy-doped Bi2212 as indicated by the triangles in figure 
10. Independent of the method chosen, the pseudogap for 
Na-doped CCOC appears significantly larger than that 
of Bi2212 at a comparable doping level. 

The high energy pseudogap seen in Bi2212 was first 
conjectured by Laughlin to be a result of the d-wave- 
like modulation of the dispersion of the insulating oxy- 
halides. 11,12 From the observed shift of the chemical po- 
tential it is clear that the high energy pseudogap indeed 
directly results from the dispersion seen in the parent 
compound, although it appears that the energy scales of 
the high energy pseudogap have a small system depen- 
dent variation. The implication for such a large pseudo- 
gap relative to the Bi2212 system at comparable doping 
should be investigated, and the bilayer effect in Bi2212 
will need to be taken into consideration. If the conden- 
sation energy for superconductivity is acquired by the 
gapping of low energy states near (n, 0), then this would 
result in superconductivity being less favorable as evi- 
denced by the smaller optimal T c (28 K and 90 K for 
Na-doped CCOC and Bi2212, respectively). 

D. Nodal Peak-Dip-Hump Structure 

A two component structure can be seen in the low en- 
ergy spectra of Na-doped CCOC at the Fermi crossing in 
the nodal direction. This was shown in figure 7e, where 
we plotted EDCs along the Fermi arc and then towards 
("7T, 0) as indicated in the cartoon inset. The reason for 
the difference in clarity of this structure between cleaves 
is unknown, but is likely due to a variation of inhomo- 
geneities and experimental conditions which make it im- 
possible to resolve the two components that appear as 
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FIG. 11. a) Temperature dependence of 12% Na-doped 
CCOC along the nodal direction. A two component struc- 
ture is observed at &f well above T c . _E 7 = 16.5 eV 

one in the (0.437T, 0.437r) spectra of figure 5. Note that 
in general, the presence of sharp low energy features is 
a tribute to a high sample quality. The loss of the two 
component structure as one approaches (w, 0) is indica- 
tive of the observed large pseudogap, which precludes the 
identification of any similar structure in the low energy 
excitations near (n, 0). 

The observation of a two component structure in these 
superconducting samples begs the question of how it may 
be related to superconductivity Figure 11 presents tem- 
perature dependence of this feature from 20 K(« T c ) up 
to 120 K. It is clear that the peak-dip-hump like struc- 
ture is observed up to 75 K, and is finally smeared out 
at 120 K. By cycling the temperature back down to 20 
K, we ensure that the broadening at high temperatures 
is not an aging effect. Since even optimally doped poly- 
crystaline samples have only a maximum T c of 28 K, it 
is clear that the two component structure observed here 
does not turn on near T c . This is in contrast with the 
more familiar two component structure in ARPES on the 
Bi2212 system which has now also been seen in YBCO 32 
and Bi2223, 33 and does turn on with T c . 34,35 It is im- 
portant to note that, aside from the different tempera- 
ture dependence, the two component structure observed 
in Na-doped CCOC occurs along the nodal direction as 
opposed to the (n, 0) region as in the other cases. Fur- 
thermore, we note that CCOC is a single layer material 
and thus the peak-dip-hump seen here can not be a con- 
sequence of bilayer splitting as has been conjectured for 
the peak-dip-hump structure observed in Bi2212. 36 

In order to gain some insight into the origin of this 
structure figure 12 presents a MDC analysis along the 
nodal direction of a x=0.10 sample. The dashed curves 
are the single Lorentzian fits to selected MDCs. Simi- 
lar procedure in other cuprates has been used to extract 
self energies from the fit parameters. 13,37,38 Here the fits 
serve as a method for parameterizing the observed dis- 
persion. Figures 12b and 12c plot the peak positions and 
widths as a function of the binding energy of the cor- 
responding MDCs. A change of slope in the dispersion 
of the peak positions is seen near 55 meV. This value 
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FIG. 12. a) Sample MDCs of Cai. 9 Nao.iCu0 2 Cl2 from the 
nodal direction taken with He la radiation. Single Lorentzian 
fits to selected MDCs are overlayed with a dotted line. Re- 
spectively, b) and c) give the peak positions and widths (half 
width at half maximum) from the MDC analysis. The arrow 
in b) indicates the kink energy, as discussed in the text. 

is attained from the intersection of two linear fits rang- 
ing from to 30 meV and 100 to 200 meV. We shall 
refer to this feature as a "kink" . A similar kink in other 
cuprate systems has received significant attention in re- 
cent ARPES literature. 13 ' 38-40 Certainly, the similarity 
of the kink energy position to the dip energy position seen 
in the EDCs suggests that these two results are related. 
The MDC width also displays an accelerated decrease in 
width as indicated by the dotted lines in figure 12c. 

What is the possible origin of a two component struc- 
ture in the EDC near Uf along the nodal direction? Pure 
macroscopic phase separation of metallic and undoped in- 
sulating domains can be ruled out, since the feature seen 
in the insulator would not have shifted to the chemical 
potential in this case. One possibility is that a sharp, 
coherent-like peak is growing with doping, while the in- 
coherent features at half filling slowly vanish. This idea 
is captured in the recent phase string calculations done 
by Muthukumar, Weng, and Sheng. 41 In this sense, the 
two components would hint at a balance between the an- 
tiferromagnetic insulator and the drive for the system 
to become metallic. An alternative scenario is that the 
peak-dip-hump structure is the result of coupling to a 
collective mode. This could also naturally explain the 
presence of the kink. Although still a debated interpre- 
tation, Lanzara et al. propose that the universally ob- 
served kink is due to electron-phonon coupling. 13 In this 
regard, it is interesting to note that the phonon breathing 
mode is expected to be roughly 10 meV less for CCOC 
on the basis of its a-axis lattice constant than in the case 
of insulating Bi2212 or LSCO as determined by optics, 42 
while the kink energy position determined by ARPES is 
less by roughly the same amount. 

Distinguishing between these two scenarios will be dif- 



ficult. A problem with the first picture would be that 
the kink would have to be explained as an artifact of 
the MDC analysis due to the multiple electronic features 
which are not individually resolved. This is not appealing 
considering the apparent universality of the kink feature 
among cuprates, 13 and one would then favor a coupling to 
a collective mode. However, it is surprising that the kink 
is seen in both LSCO where the states were created inside 
the gap, and in Na-doped CCOC where the states appear 
to originate from the effective lower Hubbard band which 
has now been shifted to the chemical potential. This sug- 
gests that a third possibility, combining both ideas, may 
also be considered for the doping evolution. Namely, that 
the chemical potential has shifted upon doping, the pre- 
cise location would be system dependent, while states 
are created at E^ which then couple to some collective 
excitations creating the universally observed kink (see 
figure Id). If this is the case, the Ca2- a; Na a ;Cu02Cl2 
and La2_ 2; Sr 2 ;Cu04 system would no longer appear so 
different. Unfortunately, the present data can not dis- 
tinguish between these three possibilities, and we must 
leave this as an opportunity for future experiments with 
more available doping levels to investigate. 

IV. CONCLUSIONS AND DISCUSSION 

The data from the valence band, near Ej? spectra, 
and photon energy dependence in the Ca2- ;r Na a ;Cu02Cl2 
system provide convincing evidence that the chemical po- 
tential shifts to the top of the valence band upon doping 
the insulator, Ca2Cu02Cl2. Surprisingly, these results 
appear to present a different evolution across the metal 
to insulator transition than in LSCO where the chemical 
potential is pinned and states are created inside the gap 
upon doping. 5 

With the comparison of CCOC we have shown that the 
high energy pseudogap is indeed a remnant property of 
the insulator. 11 Meanwhile, by examining the low lying 
excitations we found that an attempt to use a large Fermi 
surface centered at (-7T, 7t) results in a gap structure which 
deviates from the canonical d-wave functional form. In- 
stead a small Fermi pocket about (7r/2,7r/2) as naively 
expected from a rigid chemical potential shift is a better 
description of heavily underdoped Ca 2 _ x Na 2 ,Cu02Cl2. 
This is the case when we insist that the Fermi surface can 
not terminate at a point in the middle of the Brillouin 
zone. A Fermi surface arc could also explain the data if 
this were allowable in a non-Fermi liquid picture. We an- 
ticipate this region to also possess a gap structure which 
could differentiate between competing models in the un- 
derdoped regime. Finally, a two component structure 
is also observed in the EDCs along the nodal direction. 
This appears to be related to a kink seen in the MDC 
derived dispersion, which is now a universally observed 
feature of the cuprates. 13 

The presence of features seen in the insulator within 
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the low energy spectra of the Na-doped superconductor 
suggests a strong interplay between antiferromagnetism 
and superconductivity in the underdoped regime. One 
way to think of these results is to start from the over- 
doped regime where superconductivity is the first in- 
stability encountered. Upon further undcrdoping, mag- 
netism begins to play a role which opens a large pseu- 
dogap, and suppresses superconductivity. Alternatively, 
one can start with the insulator, which forms a small 
pocket as carriers are introduced into the system. Super- 
conductivity then creates a new instability that opens a 
small gap along the Fermi surface, which we are unable 
to detect given our resolution. We look forward to fu- 
ture experiments which will clarify Ca2_ £C Na a; Cu02Cl2's 
evolution to the more "conventional" behaviour of the 
overdoped regime, and its relation to the other cuprate 
superconductors. 
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